Surface width of the Solid-On-Solid models 
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The low-temperature series for the surface width of the Absolute value Solid-On-Solid model and the Discrete 
Gaussian model both on the square lattice and on the triangular lattice are generated to high orders using 
the improved finite-lattice method. The series are analyzed to give the critical points of the roughening phase 
transition for each model. 



1. INTRODUCTION 

The solid-on-solid(SOS) model is an approxi- 
mation of the interface in three-dimensional Ising 
model. It suffers from the roughening phase tran- 
sition. This transition is considered to be of 
the Kosterlitz-Thouless (KT) type@. Hasen- 
busch et alQ determined the roughening tran- 
sition points u r in high precision by Monte Carlo 
renormalization group (MCRG) matching of var- 
ious SOS models with the body-centered solid- 
on-solid(BCSOS) model, which was solved ex- 
actly. H They found that the absolute- value solid- 
on-solid(ASOS) model and the discrete gaus- 
sian(DG) model on the square lattice belong 
to the same universality class as the BCSOS 
model, with the roughening transition points 
u r = 0.19945(1) and u r = 0.2643(1) for the ASOS 
model and the DG model, respectively. 

The free energy in the SOS models is expected 
to exhibit the singularity like 

/(it) = A{u) exp [-c(it r - u)- 1/2 } + B{u). 

In fact the free energy has this type of singular- 
ity in the BCSOS model. The low-temperature 
series for the free energy of the ASOS model was 
calculated ^,|| and the behavior of the series is 
consistent with this type of singularity. This sin- 
gularity is, however, weaker than the singularity 
of the quantities concerning the interface width 
such as the average of the n-th moment of the 
surface height < h n > or the inverse of the den- 



sity gradient at the center of the interface 

M — l - 

[< 6(h + 1/2) > - < 6(h - 1/2) >] 

These quantities are expected to behave like 



< h n >,M 



(u r — u)~ 
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The renormalization group estimate was given by 
Ohta and Kawasaki M as 

1/2 for</i 2 > 
1 for < h A > . 
1/4 for M 

In the BCSOS modelQ the critical exponents 
have the same values as these RG prediction. The 
low-temperature series were calculated for < h 2 > 
and M in the ASOS model to order u 9 by Leamy 
et alQ. But, they are too short to give some def- 
inite conclusion and longer series are expected. 
Here we calculate the low-temperature series for 
< h 2 > , < h A > and M to high orders in the 
ASOS model and the DG model on the square 
lattice and also on the triangular lattice, among 
which those in the ASOS model on the square lat- 
tice are previously reported together with those 
in the 3D Ising interface. || 

2. SERIES 

The series are calculated using the improved al- 
gorithm of the finite-lattice method.) 10- 13f 
The finite-lattice can also be applies to the model 
on the triangular lattice by mapping it to the 
square lattice with not only the nearest-neighbor 



Table 1 

The low-temperature expansion coefficients for the surface height squared < ft, 2 > of the ASOS model 
and the DG model on the square lattice and on the triangular lattice. 



square lattice triangular lattice 



n 


ASOS model 


DG model 


ASOS model 


DG model 

















1 


o 


o 


o 


o 


2 


2 


2 


o 


o 


3 


8 


8 


2 


2 


4 


40 


24 


o 


o 


5 


168 


64 


12 


12 


6 


766 


172 


4 


-16 


7 


3424 


448 


78 


90 


8 


15640 


1544 





-204 


9 


71664 


6712 


552 


812 


10 


331106 


28028 


-156 


-2304 


11 


1537392 


105520 


4050 


8028 


12 


7175356 


383744 


-2476 


-25240 


13 


33633968 


1412064 


30888 


84330 


14 


158321118 


5253608 


-29004 


-273636 


15 


748150456 


19575936 


243190 


913548 


16 


3548644024 


73097908 


-305436 


-2998968 


17 


16892249464 


273313456 


1967166 


10025034 


18 


80687694550 


1016314950 


-3050432 


-33247692 


19 


386688784544 


3742138552 


16292466 


111258024 


20 


1859048866840 


13682700296 


-29581152 


-371501832 


21 


8964629577440 


49953067648 


137708278 


1246265942 


22 


43353168530402 


182594166164 


-281836944 


-4180005708 


23 


210226581434968 




1184301792 




24 






-2656332832 





Table 2 



The roughening transition points predicted by the biased inhomogeneous differential approximants of 
the low-temperature series for the surface width. The result from the MCRG matching with BCSOS 
model is also listed for comparison. 
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< h' 2 > 


0.19945(71) 


0.2636(74) 


0.3702(19) 


0.443(75) 


<h 4 > 
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0.19946(71) 


0.2643(71) 







interaction but also the next-nearest-neighbor in- 
teraction in one direction. jy| In the standard 
finite-lattice method the series expansion of the 
free energy (and its derivative with respect to the 
source term introduced to lead to the expectation 
value of various quantities) in the infinite volume 
limit can be obtained by an appropriate linear 
combination of the free energies (and their deriva- 
tives) on the finite-size lattices. The improved al- 
gorithm is effective in such a model that the spin 
variable at each site have high or infinite discrete 
degree of freedom as in the q-state Potts model 
with a large q and the SOS models. In the im- 
proved algorithm we prepare the free energies on 
finite-size lattices with the degree of freedom of 
the spin at each site restricted to be smaller than 
in the original model. Then the free energy in 
the infinite volume limit is given by an appropri- 
ate linear combination of the free energies on the 
finite-size lattices with the restricted spin degree 
of freedom as well as of the free energies with 
the full degree of freedom on smaller finite-size 
lattices. This saves us from taking into account 
unnecessarily large fluctuation of the spin vari- 
ables, which would contribute only to the higher 
order terms than we wish. The low-temperature 
expansion coefficients for < h 2 > in each model 
are listed in Table 1. The expansion parameter 
is u = exp(— 2 J/ksT), where J is the interaction 
strength in the Hamiltonian as 

J JJ2<i 3 > \ h i ~ h I fc > r ASOS model 
\ JJ2<i,j> \hi ~ hj\ 2 for DG model . 

Those for < /i 4 > and M will be seen in the 
literature. |p| 

3. ANALYSIS 

Here we assume the universality between the 
exactly solved BCSOS model and the other SOS 
models so that the critical exponents for < h 2 >, 
< /i 4 > and M would be 1/2,1 and 1/4, re- 
spectively, which is also the prediction of the 
renormalization group arguments. Then inho- 
mogeneous differential approximants of the series 
biased for the critical exponents to have the above 



values give the values of the roughening transition 
point for each model, which are listed in Table 2. 
The data listed for M is from the approximants 
of its second derivative with respect to u, which 
give more convergent result than of itself. We see 
that the obtained values of the roughening transi- 
tion point for the ASOS model and the DG model 
on the square lattice are consistent with the pre- 
cise values given by Hasenbusch et al from the 
MCRG matching to the BCSOS model. Those 
for the models on the triangular lattice are a new 
prediction, which we do not know any data to be 
compared with by now. 
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